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We establish a relation between entanglement of a many-body system and its diffractive proper- 
ties, where the link is given by structure factors. Based on these, we provide a general analytical 
construction of multi-qubit entanglement witnesses. The proposed witnesses contain two-point cor- 
relations. They could be either measured in a scattering experiment or via local measurements, 
depending on the underlying physical system. For some explicit examples of witnesses we analyze 
the properties of the states that are detected by them. We further study the robustness of these 
witnesses with respect to noise. 



Multipartite entanglement is a remarkable property of 
quantum systems that outlines fundamental discrepan- 
cies to classical physics and occurs at the interface be- 
tween quantum information and many-body physics. A 
useful tool to verify the presence of multipartite entan- 
glement are entanglement witnesses. Their experimental 
implementations via local measurements are presented, 
e.g.,inRefs. [E IE IE EL IE E| ■ Here we introduce a method 
to detect multipartite entanglement also in a scattering 
experiment, when in general no local access to the indi- 
vidual subsystems is possible. 

In this paper we present a construction for multipar- 
tite entanglement witnesses based on linear combinations 
of operators associated with structure factors. Structure 
factors consist of two-point correlations, and have been 
widely used in condensed-matter physics. They are mea- 
surable in scattering experiments, e.g. via neutron scat- 
tering in condensed-matter systems, or via light scatter- 
ing off optical lattices. For example, structure factors 
are employed in experiments 0, Hf that investigate the 
crystal structures of particular molecules that can be de- 
scribed via Hamiltionians of spin chains. One-particle 
spectral functions are investigated for various materials 
in [E IE Also, much theoretical work on determin- 
ing the structure factor and the related spectral weights 
of particular spin chain models has been done, numeri- 
cally and analytically [13, EE EI EE El • Moreover, the 
structure factor plays an important role in the physics of 
atoms in optical lattices since it is related to the visibil- 
ity of the interference pattern 17J. Impurity contribu- 



tions to the static structure factors have been predicted 
in In this paper we show how witness expectation 
values can be determined from global diffractive proper- 
ties, via the corresponding structure factors. Therefore, 



we establish a link between multipartite entanglement 
and diffractive properties of many-body systems, based 
on two-point correlations. Let us note that expectation 
values of our witnesses can be also measured locally [l9| 
in any suitable physical system where individual subsys- 
tems are accessible, e.g. entangled photons. 

We consider a multipartite quantum system consist- 
ing of N subsystems. The construction of our entan- 
glement witnesses is based on the operator S a ^(k) := 
Ei<j e^-^SfSf , where i,j denote the i-th and j-th 
spins, ri,rj their positions, and a,j3 = x,y,z. Here, Sf 
are the spin operators and k is the wave- vector transfer, 
where we consider the one dimensional situation. The 
expectation value of this operator is the static structure 
factor (apart from a different summation and normaliza- 
tion, this form appears in, e.g., Ref. EH) 



S a ?(k) 



i<j 



(i) 



which is a function of k. The dynamic formulation of the 
structure factor describes the "response" of the probed 
system [2fj. Below, we will normalize the distance be- 
tween two neighboring spins to one. The distance be- 
tween spins is either defined via the periodic structure, 
or - in a non-periodic situation, like e.g. for entangled 
photons - via the labels of the spins. So far we have con- 
sidered general spin operators. In the following we will 
focus on the case of spin 1/2, where S a corresponds to 
the Pauli operator a a . 

An entanglement witness is defined as a Hermitian 
operator W that detects the entanglement of a state 
p if it has a negative expectation value for this state, 
(W) p = Ty( P W) < while at the same time Tr(aW) > 
for all separable states a 21, 23]. We construct a class of 



2 



entanglement witnesses using the static structure factor 
operator as 



where m; 



Ti e N. Thus, for product states we 



W(k) :=t N - E(fc) 



(2) 



where 1jv is the identity operator on the 2 -dimensional 
Hilbert space and 



E(AO = i[£(fc)+E(-fc)] 



(3) 



with 



£(*) = — i— (c*s**(k) + c y syy(k) + cj"(kj) , 



1(N, 

a e JR, Id < l. 



(4) 



Here B(N, 2) is the binomial coefficient. The restriction 
to real coefficients a ensures that the witness operator 
is Hermitian, while the choice \a\ < f ensures normal- 
ization. The meaning of the parameter k depends on 
the physical system: for the detection of entanglement 
for states of multiple photons, it just fixes a sign rule 
for the two-point correlation terms of the witness and 
has no further physical meaning. For spin chains, as 
mentioned, k is the wave vector transfer in scattering 
experiments. Here the entanglement witness can be de- 
termined via a readout of values of the structure factor 
for various k. The expectation value of the witness is 
(W(k)) = 1 — (E(fc)). The crucial point is to show that 
for all product states p^ := pi <8> pi ® . . . <E> Pn we have 
(W(k)) p N > (due to convexity the same is true for sep- 
arable states that are a convex combination of projectors 
onto product states). It is helpful to remember the Bloch 
vector form of a qubit state, 



Pi 



5(1 + ,«) 2 + K) 2 +K) 2 <i. ( 5 ) 



The product state pp should then be regarded as a 
product of Bloch decompositions ([5]). In this way we 
obtain the following bound on the expectation value of 
E(fc) for product states: 

I I = 
= w^yl E ( c - cos > 

+ c y cos (kmij )(of ctJ) + c z cos (ferny ) (of crj)^ | 

= 2 J { Cx C0S ( km ij) n x n l 

+ c y cos {krtiij)n y vP y + c z cos (ferny)?!* n-M \ 



^ mYTY) £ ( |n *iki + KiKl + kiki) < 1 



have {W{k)) p N > 0. 

Which states can be detected by W(k)l For k = 
symmetric states like the Dicke states can be detected. 
These states are defined as pure states that are a superpo- 
sition of all possible permutations of I excitations (states 
|1)) in N particles and denoted as \N,l). Examples are 
the Estate \W) = |3, 1) = 1/V3(|001 ) + 1010) + 11 00)) or 
1 4, 2) = 1/V6(|0011) + |0110) + |1100) + |1001) + |1010) + 
1 0101)). Dicke states are detected by the witness ((2|) with 
c x = c y = 1, c z = — 1. To sec this, we calculate the expec- 
tation value (E(0)) of Eq. (U) for the Dicke states |JV, I). 
Entanglement is detected if (N, Z|E(0)|JV, I) > 1. For the 
term S zz {0) we get (N,l\S zz (0)\NJ) = (4(j|) -N)/2 = 
((N - 2l) 2 - N)/2. Here we have used the fact that the 
collective spin operator J a :— ~t/2J2k=i a k ^ s gi ven by 
S aa (0) = (4J2 - Nl N )/2. For S xx (0) and S yy (0) we 
obtain (^^15^(0)1^,0 = (NJ\Syy{0)\N,l) =l(N-l). 
Since for the Dicke states we have ^(O) = ^^^(O) > it 
is clear that the chances to detect entanglement are best 
for c x = c,y = \. For the case c x = c y = l,c z = — 1 we 
use the notation E(0) and find 



(N,l\t(0)\N,l) - 



H(N-l) -(N- 2l) 2 + N 
N(N - 1) 



(7) 



(6) 



In particular, for an even particle number N and I = N/2 
the expectation value becomes (E(0)) = (N+1)/(N- 
1) > 1 and for odd N and I = {N-l)/2 or I = (N + l)/2 
we get (E(0)) = (N(N + 1) - 2)/N(N - 1) > 1, and thus 
these states are always detected. Also other Dicke states 
are detected, e.g. |6,2) where (6, 2|E(0)|6, 2) = 17/15 > 
1. Choosing different coefficients c XiViZ in the construc- 
tion of the witness ([2]), other states can be detected. An 
interesting example for four particles is the superposition 
between two Greenberger-Horne-Zeilinger (GHZ) states, 
(cos0)/2(|OOll) + |1100>) ± (sin6»)/2(|0000) + |1111». 
This state is detected for 7r/4 < 8 < tt/2 if we choose 
c x = -l,c y = c z = 1 (minus sign) or c x = l,c y = 
— l,c 2 = 1 (plus sign). In Ref. [23| an experimental 
preparation of a four-qubit cluster state is reported, and 
from the presented method it seems likely that also the 
above GHZ superposition states can be prepared with 
this setup. Other detected symmetric states are super- 
positions of Dicke and GHZ states, e.g., for four parti- 
cles cos 61 1 4, 2) ± sin fl/v^f |0000) + |1111>) is detected for 
arccos3i/2/19 < 8 < tt/2 with the witness coefficients 
c x = —l,c y = c z = 1 (minus sign) and c x = l,c y = 
-l,c z = l (plus sign). 

So far we have considered the case k — only. If we 
choose k — ir in the construction of the witness W(k) 
([2]), still more entangled states can be detected. Note 
that in this case the witness is no longer symmetric un- 
der particle exchange. An example of detected states 
are non-symmetric Dicke states with additional phases. 
Choosing W(tt) with c x = c y = c z = 1 we can detect the 
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four-particle entangled state 

\Df) = i(|0011) + |1100) + |0110) + |1001) 
v6 

— 10101) - |1010)). 



(8) 



For six particles and 1^(71") with c x — c y — c z — 1 again, 
the state 

\D ph ) = 

' ( v |111000) + |001110) + |010101) + 1011010) 

+ 1 100011) + |100110) + |101001) + |101100) 
+ 1 110010) + 1001011) - |000111) - |110001) 

- 1 101010) - 1 100101) - |011100) - |011001) 

- |010110) - |010011) - |001101) - |110100)) (9) 

is detected. In general, all "phased" Dicke states \N, l ph ), 
i.e. Dicke states with different signs before terms that 
correspond to even and odd permutations of 0s and Is, 
and with I = N/2 for even N and / = (N + l)/2 or 
I = (N — l)/2 for odd N, are detected. Consider, e.g., 
the state \6,3 ph ) = \D ph ) ©. Starting with the state 
1 111000), all even permutations, i.e. an even number 
of transpositions of neighboring spins, have a positive 
sign, and all odd permutations have a negative sign. 
This scheme can be generalized to construct the general 
phased Dicke states \N,l ph ) of N particles and I excita- 
tions. 

It is not necessary that all coefficients c X:V:Z are non- 
vanishing. The states \N,l ph ) for the mentioned par- 
ticular values of I are already detected by the wit- 
ness W(it) with c x = c y = l,c z = 0. This can be 
seen as follows: To determine (N, l ph \ S xx (n)\N,l ph ) and 
(N,l ph \S yy (T:)\NJ ph ) we note that we get exactly the 
same expressions as for the Dicke states, since the possi- 
ble minus signs of the two-point correlation terms of <r* 
and the minus signs of the phased Dicke states always 
exactly cancel out, due to the correspondence of an odd 
number of transpositions in the state and odd distance- 
terms Tj — Ti in the structure factors. Thus for the phased 
Dicke states \N,l ph ) we obtain (S xx (n)} = (S yy (ir)} = 
l(N — I), just as for the Dicke states \N,l) in the pre- 
vious paragraph. For phased Dicke states with even N 
and I = N/2 and for £(7r) with c x = Cy = 1, c z = 
we get (E(tt)) = 2(S xx (tt)) = N/(N - 1) > 1, and for 
odd N with I = (N + l)/2 or I = (N - l)/2 we obtain 
(£(tt)) = (N + 1)/N > 1. Therefore these states are 
always detected. 

In order to prove that a given state carries genuine 
multipartite entanglement, one has to determine its min- 
imal expectation value for any biseparable cut of the 
multi-qubit states. This task can in general only be 
performed numerically. Using the routines provided in 
Ref. [241 ] . we find the following results for the witnesses 



for phased Dicke states, when choosing k = n and 
c x = c y = c z = 1: for 4 qubits, (Y,(tt)) bisep < 1.187 
and for 6 qubits, (£(7r))t>i S ep < 1.158, where these are 
upper bounds for all biseparable cuts. Comparing these 
numbers with the previous paragraph shows that there 
is a considerable range of genuine multipartite entangled 
states to be detected with our method. 

The importance of the parameter k for the detection 
of different types of states is remarkable: The general 
phased Dicke states \N,l ph ) that include the states (|8| 
and ((§]), are not detected by W(0) for any choice of the 
coefficients c x , c y , and c z , and W(ir) does not detect the 
"usual" Dicke states. 

Furthermore, we want to study the robustness of the 
witness in Eq. ^ under the influence of noise. In partic- 
ular, we consider two depolarizing channels; one that acts 
collectively on all qubits, which efficiently results in the 
addition of white noise to the multipartite state, and one 
that affects the single qubits independently, i.e. that adds 
white noise to a single qubit p s ; p s ,dp = (1 — <z)Ps + <7l/2- 

We study the collective depolarizing channel first. 
The Dicke states \N,N/2) then change accordingly to 
pl N /2 N + (l-p) \N,N/2)(N,N/2\. The witness W(0) 
with c x — c y — 1, c z = — 1 detects entanglement of this 
state for < p < 2/(N +1). The robustness decreases 
with a growing number of qubits N. In Ref. 25j, the 



robustness of certain witnesses against the collective de- 
polarizing channel, corresponding to c z = 0, has been 
studied. These witnesses, that can detect Dicke states 
of the form \N,N/2), allow noise with p < 1/N. Thus, 
adding the z-direction measurement to the witness im- 
proves its robustness against white noise. For the general 
"phased" Dicke states and using the witness W(ir) with 
c x = c y = 1 and c z = 0, we obtain < p < 1/N for 
entanglement detection. Again, a greater robustness can 
be achieved when an additional z-direction term, c z — 1, 
is introduced. For the noisy "phased" Dicke state (cf. 
Eq. (JU)), i.e. pl/16+(l-p)\D Ph ){D Ph \, the witness W{n) 
with c x = c y = c z = 1 detects entanglement of this state 
for < p < 4/13. In the case of six particles, it detects 
entanglement of the noisy state pi/ 2 6 + (l—p)\D Q h )(Dg h \ 
(cf. Eq. ((5J)) for a parameter interval < p < 6/31. The 
maximal values of p are in both cases bigger than 1/N. 

In the following we investigate the robustness for the 
individual depolarizing channel, where the noise affects 
each qubit independently. Since we are interested in ex- 
pectation values only, it is convenient to shift the influ- 
ence of the Kraus operators K+ characterizing the noise 
model to the observable and leave the initial state un- 
changed. This is possible, because in the operator sum 
representation of the channel we have (where the sub- 
script dp denotes affection by the channel) Tr (Opd p ) = 
^TiiOKipKj) = ^TiiKjOKip) = Tr(O dpP ). Since 
the individual depolarizing channel transforms the Pauli 
operators as cr^ = (1 — q)a a , a — x, y, z (see, e.g., [26|), 
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the two-point correlation terms of the structure factor 
simply change to (o-fo^)dp = (1 — q) 2 o~i • Thus the 
observable S(fc) in Eq. (H|) is influenced by the channel 
according to T,(k)d p = (1 — q) 2 £(fc), and for the expec- 
tation value we simply have (£(fc)<£ p ) = (1 — q) 2 (I](k)). 
Thus we can determine the robustness of W(k), i.e. the 
region of q for which (^(fc)^) > 1, where entangle- 
ment is detected. For the Dicke states \N, N/2) and 
c x = c y = 1, c z = — 1 for W^(0) we obtain a robustness 
region < q < 1 - \f(N — 1)/(N + 1). For the phased 
Dicke states \N, N/2P h ) and c x = c y = l,c z =0 for W(tt) 
we find < q < 1 - y/(N-l)/N. For c x = c y = c z = 1 
we calculated the four and six particle case (see Eqs. ((HI) 
and ©), where we find < q < 1 - 3/VT3 ~ 0.168 for 



\D 



ph v 



and < g < 1 - 5/V31 ~ 0.102 for \D, 



iP h \ 



Both 



states exhibit a robustness region that is larger for c z = 1 
than for the c z = case, it therefore seems favorable to 
add the z-direction measurement setting also in the case 
of the individual depolarizing channel. 

While writing this paper, interesting experimental and 
theoretical articles on detecting multi-qubit entangle- 
ment of photonic Dicke states appeared online [1, H, [27], 
HU . The witnesses that were discussed and implemented 
there can be seen as specific cases of the general con- 
struction presented here. A detailed noise study for wit- 
nesses of Dicke states that includes the amplitude and 
phase damping channels, in connection with the numer- 
ical detection of genuine multipartite entanglement, can 
be found in Ref. [281 ] . Note that a different method for 
macroscopic detection of non-separability, which is valid 
for a certain class of spin models, is provided via the 
magnetic susceptibility, as shown in Ref. [29j ], 

Summarizing, we present a general construction of en- 
tanglement witnesses for multi-qubit states and establish 
a relation to static structure factors - these are macro- 
scopic quantities which can be determined or measured 
in a collective way in various periodic physical systems 
such as spin chains. Our approach opens a wide avenue 
of possible new connections for experimental detection 
of entanglement through scattering, and can lead to a 
deeper understanding of entanglement properties in con- 
densed matter systems. Explicit constructions of wit- 
nesses for higher spins and eventually continuous variable 
systems can follow along the same lines. A connection to 
dynamical structure factors may allow also to study the 
dynamics of entanglement. Our method can be tested on 
spin chains with a finite number of sites, as for example 
in molecular magnetic materials (30j . 

In suitable physical systems where the constituents 
can be addressed individually (such as polarized pho- 
tons) the witnesses can be measured locally. Our wit- 
nesses are suitable for detecting a diversity of states, such 
as Dicke states, non-symmetric versions of Dicke states, 
GHZ states and superpositions of these. The presented 
construction also allows variations: in principle, arbitrary 



signs of the constituting two-point correlations are pos- 
sible. This could be useful for the detection of various 
entangled states in photon experiments. Thus, we offer 
more generality than previous two-point entanglement 
witnesses like spin squeezing inequalities (see Ref. 31 1 
for an overview). The question of whether the witnesses 
(that are constructed for a finite number of subsystems) 
are also meaningful in the thcrmodynamical limit is still 
under investigation. 
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